2, 2008 3:31 WSPC/INSTRUCTION FILE kniehl 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



HEAVY-QUARKONIUM PRODUCTION AT NEXT-TO-LEADING 

ORDER 

BERND A. KNIEHL 

//. Institut fiir Theoretische Physik, Universitdt Hamburg, Lumper Chaussee 149, 
22761 Hamburg, Germany 

We review recent progress in the description of heavy-quarkonium production in 2 ^ 2 
processes at next-to-leading order in the factorization framework of nonrelativistic quan- 
tum chromodynamics. Specifically, we consider the production of prompt charmonium in 
association with a hadron jet or a prompt photon in two-photon collisions and exclusive 
double-charmonium production in e+e" annihilation. 
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1. Introduction 

The factorization formalism of nonrelativistic quantum chromodynamics (NRQCD) 
provides a rigorous theoretical framework for the description of heavy-quarkonium 
production and decay that is renormalizable and predictiveP Theoretical predic- 
tions are decomposed into sums over products of short-distance coefficients, which 
can be calculated perturbatively as expansions in the strong-coupling constant as, 
and long-distance matrix elements (MEs), which are subject to relative- velocity (v) 
scaling rules (see Table ^ and must be extracted from experiment, and are so orga- 
nized as double expansions in as and v. This formalism takes into account the com- 
plete structure of the QQ Fock space, which is spanned by the states n = "^s+ijj."-) 
with definite spin S, orbital angular momentum L, total angular momentum J, 
and color multiplicity a = 1,8, and so predicts the existence of color-octet (CO) 
processes in nature. In the case of S'-wave quarkonia, the traditional color-singlet 
(CS) model^ is recovered in the limit w — > 0. However, the latter suffers from severe 
conceptual problems indicating that it is incomplete. These include the presence of 
infrared (IR) singularities in the P-wave case and the lack of a general argument 
for its validity in higher orders of perturbation theory. 

The greatest triumph of the NRQCD factorization formalism was its ability to 
correctly describe the cross section of inclusive charmonium hadroproduction at the 
Tevatron, which exceeds the prediction in the CS model by more than one order of 
magnitude. In order to convincingly establish the phenomenological significance of 
the CO processes, it is indispensable to identify them in other kinds of high-energy 
experiments as well. The verification of the NRQCD factorization hypothesis is 
presently hampered both from the theoretical and experimental sides. On the one 
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Table 1. Values of k in the veloc- 
ity-scaling rule (O^fn]) oc for the 
leading cc Fock states n pertinent to 
H = J/i>,Xcj,'>P'- 

k J/tp, V' XcJ 

3 35(1) 

5 - 3p(i)^ 35(8) 

7 15^®', ^S'®', 3pj8) _ 

hand, the theoretical predictions for direct heavy-quarkonium production to be com- 
pared with existing experimental data are, apart from very few exceptions ,I2EEEI 
of lowest order (LO) and thus suffer from considerable uncertainties. The measure- 
ment of charmonium polarization at the Tevatron currently presents a challenge 
for NRQCD factorization, but any conclusions are premature in the absence of a 
full-fledged next-to-leading-order (NLO) analysis. It is, therefore, mandatory to cal- 
culate the NLO corrections to the hard-scattering cross sections and to include the 
composite operators that are suppressed by higher powers in v. Apart from the 
usual reduction of the renormalization and factorization scale dependences, sizeable 
effects, e.g. due to the opening of new partonic production channels, are expected 
at NLO. On the other hand, the experimental errors are still rather sizeable. The 
latter are being significantly reduced by HERA II and run II at the Tevatron, and 
will be dramatically more so by the LHC and hopefully a future e'^e~ linear col- 
lider such as the TeV-Energy Superconducting Linear Accelerator (TESLA) , which 
is presently being designed and planned at DESY. 

Recently, 2 — s- 2 processes of heavy-quarkonium production were for the first 
time studied at NLO in the NRQCD factorization formalism. Specifically, the pro- 
duction of prompt charmonium, which is produced either directly or through the 
decay of heavier charmonia, with finite transverse momentum (pt) in association 
with a hadron jeJ^or a prompt photon^ via direct photoproduction in two-photon 

. . 4- _ . . . ft 

collisions and exclusive double-charmonium production in e e annihilatiorii^ were 
considered. In this presentation, we review the most important conceptional issues 
and phenomenological results of Refs. 131121 in Sections [3 and |31 respectively, report 
on Ref. iniin Sectional and offer an outlook in Section |21 

2. Conceptional issues 

We focus attention on the process 77 — s- J/iJj + X, where J/tp is promptly produced 
at finite value of px and X is a purely hadronic remainder. Since the incoming pho- 
tons can interact either directly with the quarks participating in the hard-scattering 
process (direct photoproduction) or via their quark and gluon content (resolved pho- 
toproduction), this process receives contributions from the direct, single- resolved, 
and double-resolved channels, which are formally of the same order in the pertur- 
bative expansion. At LO, the bulk of the cross section is due to single-resolved 
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photoproduction, and the NRQCD predictioiiZi based on the MEs determined from 
fits to Tevatron dat nicely agrees with a recent measurement by the DELPHI 
Collaboration at LEP2P 




Fig. 1. Tree-level Feynman diagrams for process ITl . 



Here, we consider direct photoproduction at NLo!^ At LO, there is only one 
partonic subproccss, namely 



7 + 7- 



9- 



1 



The relevant Feynman diagrams are depicted in Fig.Q] At NLO, virtual corrections 
to process and {O^^sf'^]), where H = J/ip^XcJiip'^ a-nd real corrections to 



7 + 7 — > cc[n\ + g + g, 
7 + 7 ^ cc[n] + u + M, 
7 + 7 ^ cc[n] + q + q, 



n 
n 



3p(l) lc.(8) 3q{S) 3p(8) 



■ipy-) iq\°) -icy^ 
jr J , Oq , Ji 

3o(8) 
•^1 ; 



J ' 



„ _ lq{S) 3c.(8) 3p(8) 



(2) 
(3) 

(4) 



where u and u denote the Faddeev-Popov ghosts of the gluon, contribute. 

The virtual corrections to process ^ receive contributions from self-energy, 
triangle, box, and pentagon diagrams. The self-energy and triangle diagrams are in 
general ultraviolet (UV) divergent; the triangle, box, and pentagon diagrams are in 
general IR divergent; and the pentagon diagrams without three- gluon vertex also 
contain Coulomb singularities. As for the light-quark loops, the triangle diagrams 
vanish by Furry's theorem, while the box diagrams form a finite subset. The virtual 
corrections to {0^[^s[^'']) also produce UV, IR, and Coulomb divergences. The 
UV and IR divergences are extracted using dimensional regularization in d 4 — 
2e space-time dimensions, leading to poles in eyv and eir,, respectively, while the 
Coulomb singularities are regularized by a small relative velocity v between the c and 
c quarks. The UV divergences are removed by the renormalization of (0^[^iS'|;*'']), 
as, the charm-quark mass and field, and the gluon field, which is performed in the 
modified minimal-subtraction (MS) scheme for the former two quantities, rendering 
them dependent on the renormalization scales A and fj,, respectively, and in the 
on-mass-shell scheme for the residual three quantities. The IR divergences cancel 
among the virtual and real corrections, the wave-function renormalizations, and 
(O^pS-f^]). The Coulomb diver gences cancel between the virtual corrections and 
(0^p5f)]). 
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The real corrections are plagued by IR divergences, which come as collinear 
divergences from the initial state and collinear and/or soft ones from the final state. 
They are identified by appropriately slicing the three-particle phase space using 
small parameters Si and Sf, respectively. The collinear and/or soft regions of phase 
space are integrated over analytically in d dimensions, while the hard region is 
integrated over numerically in four dimensions. The sum of these contributions is, 
to very good approximation, independent of 6i and 6f. The initial-state collinear 
divergences are factorized at some factorization scale M and absorbed into the 
parton density functions (PDFs) of the q and q quarks inside the resolved photon. 
The M dependence thus introduced is approximately compensated by the LO single- 
resolved contribution. 

Combining the contributions arising from the virtual corrections (vi), the pa- 
rameter and wave- function renormalization (ct), the operator redefinition (op), the 
initial-state (is) and final-state (fs) collinear configurations, the soft-gluon radiation 
(so), and the hard-parton emission (ha) as 

da{n, A, M) = (ifTo(A*i ^)[1 + '^vi(M; ^uv, eiR, + '^ct(M; euv, em) + Sop{^J-, A; eir, v) 
+ Sfsin; eiR, 6f)] + da-.silJ., A, M; 5j) -I- daso{n, em, Sf) + dcrha(M: ^i, ^f), (5) 

the regulators euv, eiR, u. Si, and Sf drop out and the fi and A dependences formally 
cancel up to terms beyond NLO, while the M dependence is unscreened at NLO. 

3. Phenomenological results 

We consider two-photon collisions at TESLA operating at center-of-mass en- 
ergy -y/g = 500 GeV, where the photons arise from electromagnetic initial-state 
bremsstrahlung, with antitagging angle 0max = 25 mrad, and beamstrahlung, with 
effective beamstrahlung parameter T = 0.053. The J/ip, XcJ, and ■0' MEs are 
adopted from Ref. ^ and the photon PDFs from Ref. j^l 

In Fig. 121 we study d^a/dpTdy (a) for rapidity y = as a function of pT and (b) 
for pt — 5 GeV as a function of y, comparing the LO (dashed lines) and NLO (solid 
lines) results of direct photoproduction with the LO result of single-resolved photo- 
production (dotted lines). From Fig.|2Ia), we observe that, with increasing value of 
Pt, the NLO result of direct photoproduction falls of considerably more slowly than 
the LO one. In fact, the QCD correction (K) factor, defined as the NLO to LO ratio, 
rapidly increases with px, exceeding 10 for pT > 10 GeV. This feature may be under- 
stood by observing that so-called fragmentation-pron^^partoiiic subprocesses start 
to contribute to direct photoproduction at NLO, while they are absent at LO. Such 
subprocesses contain a gluon with small virtuality, q^ = 4m^, that splits into a cc 
pair in the Fock state n = and thus generally generate dominant contributions 
at Pt ^ 2mc due to the presence of a large gluon propagator. In single-resolved 
photoproduction, a fragmentation-prone partonic subprocess already contributes at 
LO. This explains why the solid and dotted curves in Fig.|2Ka) run parallel in the 
upper Pt range. At low values oipT, the fragmentation-prone partonic subprocesses 
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Pt IGeV] y 

Fig. 2. LO single-resolved (dotted lines), LO direct (dashed lines), and NLO direct (solid lines) 
contributions to d? cr / dpq^dy (a) for i^ = as a function of pj- and (b) for = 5 GeV as a function 
of y. 



do not matter, and the relative suppression of direct photoproduction is due to the 
fact that, at LO, this is a pure CO process. 

4. Exclusive double-charmonium production in e+e~ annihilation 

Until recently, another serious challenge for NRQCD factorization emanated from 
the fact that the exclusive cross section of e+e^ J/tp + rjc measured by the 
BELLE Collaboratioiff2l at KEKB with ^ « 10.6 GeV significantly exceeded the 
LO predictionpl by a factor of 5-10. Fortunately, the situation was mitigated by a 
NLO analysis^ and an independent measurement by the BABAR Collaboration at 
PEP-II.EI 

Through NLO, e+e" J/tp+rjc is a pure CS process.l^lThe NLO corrections are 
entirely virtual and receive contributions from self-energy, triangle, box, and pen- 
tagon diagrams. The UV divergences are removed by parameter and wave-function 
renormalizations, the IR ones cancel against those from the wave-function renor- 
malizations, and the Coulomb ones are absorbed into the CS MEs (O'^/'^psf ^]) and 

In Table 21 the NLO result for the total cross section ^(e+e — > J/ip -\- rjc) 
at = 10.6 GeV including the theoretical uncertainty from t he c hoice of renor- 
malization scale is compared with the BELLEiEl and BABAbSSI 

measurements, 

where the first errors are statistical and the second ones systematic. Notice that the 
experimental values include the branching fraction of the rjc decays into final states 
containing more than two charged tracks. 



5. Outlook 

In order to complete the NLO treatment of prompt J / tjj production in two-photon 
collisions, one still needs to evaluate the NLO corrections to single- and double- 
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Table 2. a{e^e —> J/ip + rjc) (in fb) as pre- 
dicted in and measured by BELLE*^^ and 

babarEI 



NLO 



BELLE 



BABAR 



14.1 



+8.9 
-3.8 



> 25.6 ± 2.8 ± 3.4 



> 17.6 ± 2.8 



+ 1.5 
-2.1 



resolved photoproduction. Then, also prompt J/tp production in photoproduction 
at HERA and hadroproduction at the Tevatron can be described at NLO. Further- 
more, in order to fully exploit the measurements of exclusive double-charmonium 
production at the B factories, also other pairs of charmonium states should be con- 
sidered at NLO. These studies will provide a solid basis for an ultimate test of the 
NRQCD factorization framework. 
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